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ABSTRACT 


Angular  scattering  diagrams  are  constructed  for  spherical  aerosols  of  refrac¬ 
tive  indox  n  =  1.  5.  They  show  the  intensity  functions  ip  iy,  and  i^  +  i^  as 
functions  of  the  scattering  angles  for  size  parameters  0.  5  (0.  5)10.  The 
basic  data  have  been  computed  using  Mie's  theory,  and  a  graphical  interpola¬ 
tion  technique  (altitude -chart  technique)  has  been  used  to  determine  accurately 
the  position  and  numerical  value  of  the  maxima  and  minima. 

The  scattering  diagrams  presented  in  this  atlas  allow  to  determine  the  intensity 
functions  ip  i2»  *1  +  ^2  scattering  angle  with  a  high  degree 

of  reliability.  The  basic  data  limit  the  construction  of  reliable  scattering 
diagrams  to  a<  10. 
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I.  INTRODUCTION 


Theoretical  data  for  angular  acattering  functions  for  spherical  aerosols  are  needed 
for  refractive  indices  n  between  1.0  and  2.0.  For  n  =  1.33,  scattering  diagrams 
have  been  published  in  an  earlier  report  (Penndorf,  1961b).  Herein,  samples  are 
given  in  the  form  of  an  atlas  for  n  =  1,5,  so  that  the  user  can  take  the  needed  in¬ 
formation,  from  the  curves,  which  is  more  reliable  than  a  numerical  interpola¬ 
tion  of  tabulated  values.  The  interpolation  technique  is  the  same  as  in  the  last 
report.  It  had  occurred  to  us  that  the  samples  given  in  steps  of  Aa  =  0.5  may  not 
be  sufficient:  hence,  a  large  number  of  workers  in  this  field  have  been  asked  as 
to  their  preference  of  the  steps  to  be  selected.  It  turns  out  that  for  most  practi¬ 
cal  application  an  atlas  in  steps  of  Aa  =  0.5  seems  to  be  sufficient;  however,  for 
more  detailed  evaluations,  the  original  tables  are  desirable.  It  is  planned  to 
publish  these  in  the  future  also. 

Computations  of  scattering  functions  for  isotropic  spherical  nonabsorbing  aerosols 
are  based  on  Mie's  theory,  which  leads  to  exact  numerical  values.  The  basic 
quantities  are  the  intensity  function  i^  and  i2.  They  have  the  following  meaning: 
the  intensity  function  ii  and  i2  are  proportional  to  the  magnitude  of  the  two  in¬ 
coherent,  plane-polarized  components  scattered  by  one  particle;  i^  is  propor¬ 
tional  to  the  intensity  whose  electric  oscillations  are  perpendicular  to  the  plane 
of  scattering:  i.e. ,  polarized  light  in  the  horizontal  plane:  i2  is  proportional  to 
the  intensity  whose  electric  oscillations  are  in  the  plane  of  scattering:  i.e. , 
polarized  light  in  the  vertical  plane. 

The  number  of  diagrams  in  the  atlas  is  limited  to  the  range  a  =  0.5  to  10.0,  be¬ 
cause  the  reliability  for  interpolation  decreases  above  a  =  10,  as  has  been  dis¬ 
cussed  in  the  last  report,  so  that  it  does  not  seem  justified  to  construct  diagrams 
for  large  a  values  although  numerical  data  exist  up  to  a  =  30  for  scattering  angles 
1  (1)  10  (10)  180*  (Penndorf-Goldberg,  1953). 
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11.  DATA 


The  atlae  ie  baaed  on  numerical  valuea  for  the  intensity  functions  ii  and  which 
have  been  computed  using  an  IBM  701  (Penndorf-Coldberg.  1953).  They  are 
available  for  a  refractive  index  n  =  1.5,  size  parameters  a  s  0. 1  (0. 1)  30,  and 
scattering  angles  6=0(1)  10  (10)  180*,  where  6=0*  defines  the  scattering  angle 
in  the  forward  direction. 

A  detailed  table  containing  ij  and  i2  in  the  same  steps  of  A6  but  for  a=  0.2  (0.2) 
159  has  been  published  by  Giese  et  (1962)  and  used  for  comparison  and  detail 
above  a  =  30  in  the  range  6  =  0  to  10*.  The  listed  values  for  ii  and  i2  agree  very 
well  with  our  table.  In  addition,  Olaf  and  Robock  (1961)  and  Robock  (1962)  have 
published  scattering  diagrams  for  n  =  1.5  and  a=  2,  8,  10,  20,  26,  and  30.  Their 
basic  computations  have  been  quite  detailed;  namely,  steps  of  A6  =  2*  for  oe  <  10 
and  A6  =  1*  for  a  >  12.  The  scattering  diagram  for  ot  =  8  has  been  found  useful, 
that  for  a  =  10  is  very  difficult  to  read  and  could  be  used  only  partially. 

The  investigation  of  our  data  shows  that  the  steps  of  Aoe  =  0. 1  in  the  size  param¬ 
eter  are  sufficient,  but  the  choice  of  scattering  angles;  namely,  a  constant  step 
of  A6  s  10*,  was  unfortxmate.  While  this  step  is  sufficient  to  describe  the  in¬ 
tensity  functions  up  to  about  a  =  3,  it  is  in  general  insufficient  for  larger  values 
of  ot  as  far  as  simple  numerical  interpolation  is  concerned.  Only  with  the  help 
of  sophisticated  graphical  interpolation  techniques  is  it  possible  to  extend  the 
useful  limit  of  the  tabulated  values  to  about  a  =  10.  It  is  not  recommended  to 
interpolate  our  tables  for  a  >  10  except  for  the  first  diffraction  minimum. 
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III.  RESULTS 


3.  1  ALTITUDE  CHART  FOR  INTENSITY  FUNCTION  ij 

The  interpolation  methods  have  been  discussed  in  detail  in  the  last  report. 

Hence,  there  seems  to  be  no  need  for  repeating  them.  The  altitude  charts  have 
been  constructed  for  i^  and  and  the  basic  features  appear  similar  to  those 
for  n  =  1.  33.  The  altitude  chart  for  i,  for  the  range  of  a  <  10  is  shown  in  fig¬ 
ure  3-1,  for  the  range  20  <  a  < 40  in  figure  3-2.  For  the  range  10  <  a  <  40,  an 
altitude  chart  for  0*  <  0  <10"  has  been  plotted. 

Two  completely  different  types  of  contours  appear  on  figure  3-1,  one  for  0  <90* 
and  another  for  90  <0  <  180*.  In  the  forward  scattering  area,  large,  wide,  and 
flat  ridges  are  found,  the  peaks  of  which  (indicated  by  the  letter  H)  occur  at 
larger  p  values*  the  larger  u.  For  u  =  0,  a  peak  occurs  at  p  =  6;  for  u  =  3.  6,  at 
p  =  6.  5  and  7.  5;  and  for  u  =  6.  0,  at  p  =  8.  8;  while  the  fourth  peak  is  situated  out¬ 
side  of  this  chart.  The  peaks  are  situated  on  ridges,  which  are  cut  off  by  quasi¬ 
vertical  running  valleys  into  which  deep  trenches  are  embedded.  The  valleys 
show  some  meandering  around  the  classical  diffraction  values;  namely, 
u  =  3.  83,  7.  02,  and  10.  17,  which  they  approach  asymptotically  for  very  large 
values  of  p.  Deep  trenches  occur  systematically  in  each  valley.  Investigating 
the  first  diffraction  minimum  up  to  a=  40  (figures  3-1  and  3-2),  which  can  be 
done  because  data  for  i^  and  ^2  are  available  in  steps  of  Ad  =  l*up  to  0  =  10*,  it 
was  found  that  these  trenches  occur  at  p  =  7.  2,  10.  6,  19.6,  23.4,  26.  3,  29.8, 

32.  2,  36.  1,  and  38.  8.  The  range  10  <  a  <18  is  difficult  to  evaluate,  and  the 
trench  position  is  not  very  accurate.  Therefore, we  have  not  estimated  any  values 
in  this  range  now  but  will  do  so  as  the  next  step.  Hence,  the  trenches  are 
formed  at  intervals  of  Ap=sJt»  whenever  u  reaches  an  extreme  value  (maximum 
or  minimum).  This  again  is  an  important  result  because  the  position  of  the 
deepest  points  of  the  trenches  can  be  estimated  without  additional  computations. 

The  position  of  the  trenches  has  been  investigated  further.  For  p  >  6,  the  first 
diffraction  minimum  occurs  at  scattering  angles  0  <  20  .  The  cross  section 
theorem  relates  the  total  scattering  coefficient  K  and  |o|  >  Hence,  the  question 
arises  whether  a  useful  relationship  exists  between  K  and  ij  |l®^  min(  .  It  turns 
out  that  K  reaches  a  maximum  or  minimum  at  the  same  size  parameter  pat  which 
the  deepest  point  occurs  in  the  trenches  of  the  first  diffraction  minimum.  This 
particular  p  value  is  independent  of  n  in  the  range  1.  0  <  n  <  1.  5.  It  has  not  been 
checked  for  larger  n  values,  but  it  is  believed  that  this  rule  is  valid  up  to  about 
n  =  2. 


*pis  the  normalized  size  parameter  and  defined  as  p-  2a(n  -  1).  In  the  case 
of  n  =  1.  5,  p=a.  The  diffraction  parameter  is  defined  in  the  last  report. 
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The  results  are  shown  in  table  3*1.  The  running  integer  y  indicates  the  first, 
second. . .  trench.  The  first  group  shows  the  trench  position  related  to  u<  3.  83 
when  K  reaches  a  minimum,  the  second  group  contains  those  where  u  >  3.  83 
and  K  reaches  a  maximum.  The  position  for  K  when  ne,l.  0  is  given  in  the  first 
column  (after  Penndorf,  1958),  the  second  column  the  trench  position  for  n  =  1.  33 
(Penndorf,  1961a),  and  the  last  column  for  n  =  1.  5  as  obtained  from  figures  3-1 
and  3-2  and  unpublished  parts  of  our  altitude  charts.  Itis  seen  that  the  position 
of  the  trench  is  systematic  and  can  be  expressed  as 

P  j  =  2ff(y  +.  1/4), 

where  pj  is  the  position  of  a  trench  (or  a  maximum  and  minimum  in  K).  The 
plus  sign  refers  to  u<  3.  83  and  a  minimum  in  K;  the  minus  sign  refers  to 
u  >  3.  83  and  a  maximum  InK.  Interpolated  values  are  shown  in  parentheses  (  ). 

In  the  case  when  K  reaches  a  maximum, the  agreement  is  very  good;  in  the  case 
of  K  reaching  a  minimum,  the  trench  occurs  a  bit  earlier  (ZipesO.  3-0.  5).  The 
agreement  for  y  =  1  is  always  poor,  as  expected. 


TABLE  3-1 

POSITION  OF  TRENCHES  FOR  FIRST  DIFFRACTION  MINIMUM 


K  . 
min 

and  u  <3.83 

^max  >  3.  83 

K 

K 

n 

4 

y 

n  =  1.  0 

1.  33 

1.  5 

n  =  1.  0 

1.  5 

1 

7.  6 

6.  9 

7.  2 

4.  1 

— 

— 

2 

14.  0 

13.  5 

(13.  5) 

10.  8 

9.  9 

10.  6 

3 

20.  3 

19.  8 

19.  6 

17.  2 

16.  8 

(17.  0) 

4 

26.  6 

26.  3 

23.  5 

23.  4 

5 

32.  9 

32.  2 

29.  8 

29.  8 

6 

39.  2 

38.  8 

36.  1 

36.  1 

Pt 

2ff(y  +  1/4) 

The  table  lists  the  value  pat  which  the  maxima  and  minima  in  K  occur  for  u  =  1.  0 
and  the  deepest  point  of  trenches  forij  and  n  =  1.  33  as  well  as  n  =  1.5, 
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For  the  second  diffraction  minimum,  the  first  trench  is  found  at  p=  6.  5  and 
u  =  5,  5,  and  the  next  at  p=  9.  0.  Because  of  the  interval  =  10*  in  this  domain 
of  the  altitude  chart,  the  last  one  cannot  be  exactly  located. 

Other  diffraction  minima  have  also  been  identified  and  can  be  seen  on  the  chart. 
Basically,  the  system  is  similar  to  that  of  n  =  1.  33  and  the  minima  start  between 
6  =  80  and  90*and  move  toward  smaller  angles  as  p increases. 

However,  there  are  differences  in  the  charts  for  n  =  1.  33  and  1.  50,  which  are 
based  on  a  more  careful  investigation  for  n  =  1.  40,  1 . 44,and  1. 50,  and  which  have 
led  to  several  improvements.  One  point  is  the  merging  of  diffraction  minima. 
Whenever  the  first  diffraction  reaches  a  trench  where  u  is  a  maximum,  the  first 
and  second  diffraction  minima  merge.  This  is  seen  on  figure  3-1  at  u  =  4.  6  and 
p=  9.  5,  and  even  better  on  figure  3-2.  In  the  area  p>  23,  both  diffraction 
minima  come  close  together  and  run  parallel  for  a  while.  The  two  valleys  are 
definitely  separated  by  a  ridge  of  higher  values  with  a  peak  around  p=  24.  6.  It 
is  drawn  as  if  the  first  minimum  ends  near  p=  25.  5,  but  it  could  have  been  con¬ 
tinued  and  both  connected  at  p=  26.  3.  From  a  practical  point  of  view,  it  is  un¬ 
important  whether  they  actually  merge  or  run  parallel  for  a  while.  The  distance 
between  the  two  minima  is  about  1*  for  p=  20  to  30  and  could  not  be  distinguished 
on  a  scattering  diagram.  The  same  thing  happens  again  around  p=  29.  8,  where 
both  minima  run  parallel  for  a  while  finally  to  merge  at  about  p=  31.  4.  From 
this  detailed  study,  based  on  figure  3-2,  it  is  concluded  that  the  actual  merger 
of  minima  takes  place  and  our  interpretation  at  p=  9.  5  is  probably  correct.  The 
consequence  of  such  mergers  is  a  gradual  reduction  of  the  number  of  minima 
indicating  the  incorrectness  of  our  earlier  extrapolation  (table  I,  Penndorf,  1961b). 
Rowell  (private  communication)  has  completed  studies  for  n  =  1.  486  in  steps  of 
l*and  his  interpretation  does  not  show  mergers  as  we  have  discussed  them, 
but  a  termination  of  the  first  minimum  whenever  it  reaches  its  lowest  value  of  u. 

It  seems  important  to  us  to  check  this  point  further  and  clarify  the  interpretation 
of  the  data  to  see  whether  our  interpretation  or  that  of  Rowell  is  correct. 

The  backscatter  area  (&>90°)  shows  an  entirely  different  pattern  of  the  contours. 
Deep,  long  valleys  occur  in  nearly  straight  lines  into  which  deep  trenches  are 
embedded  at  regular  intervals.  The  valleys  meander  somewhat  around  a  straight 
line,  but  much  less  than  the  diffraction  minima.  The  position  and  depth  of  the 
valleys  are  more  difficult  to  establish  than  for  n  =  1.  33,  because  the  selected 
spacing  of  Afl  =  lO'is  too  wide  for  large  a.  However,  basically,  the  interpolation 
problem  up  to  q=  10  is  not  too  difficult.  The  distance  between  the  centers  of  the 
trenches  for  a  particular  valley  is  about  8.  This  value  is  somewhat  larger 

than  for  n  =  1.  33,  where  it  was  found  to  be  about  0.  65  to  0.  75. 

For  p=  8.  0,  our  data  can  be  compared  with  the  scattering  diagram  published  by 
Olaf  and  Robock  (1961),  and  it  is  believed  that  they  agree  more  or  less  with  our 
evaluation  (figure  3-1  and  the  appendixes).  However,  the  curve  for  a=  10.  0 
published  by  Robock  (1962)  is  difficult  to  evaluate  and  could  not  be  used. 
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The  minima  shown  in  figure  3-1  are  redrawn  on  figure  3-3a,  where  the  dots 
indicate  the  position  of  the  trenches.  It  seems  possible  to  us  to  explain  the 
position  of  trenches  as  well  as  the  system  of  the  minima  (dark  rings).  In  our 
last  report,  we  have  already  indicated  a  plausible  interpretation  of  the  position 
of  the  dark  rings  for  i,  &nd  n  =  1.  33  as  a  combination  of  diffraction  minima  and 
reflection  minima.  This  line  of  investigation  is  now  pursued  further.  The  term 
reflected  ray  is  used  for  those  rays  which  hit  the  spherical  particles  and  undergo 
internal  reflections  and  refractions  at  the  surface  (see  figure  37  in  van  de  Hulst, 
p.  201  (1957)  ). 

A  very  detailed  and  careful  analysis  of  the  original  data  in  the  altitude  chart 
(figure  3-1  only  shows  selected  isolines)  leads  to  two  systems,  which  are  in¬ 
dicated  on  figure  3-3b.  The  solid  line  system  (diffraction  minima)  always  starts 
at  8  =  180°,  it  moves  toward  the  90*  area  more  or  less  at  a  constant  p  value, 
and  then  with  a  sharp  knee  downward.  The  second  system  of  dashed  lines  (re¬ 
flection  minima)  runs  more  or  less  parallel  to  the  6=  180* line.  These  two  systems 
intersect  each  other  at  the  trenches  (dots  on  figure  3-3b,  which  are  identical  with 
those  on  figure  3-3a).  We  interpret  this  finding  in  the  following  way.  The  solid 
lines  represent  the  minima  caused  by  diffraction  and  the  dashed  lines  those  by 
internal  reflection  and  refraction.  Whenever  these  two  systems  cross  each 
other,  a  trench  is  created.  This  seems  plausible  because  the  trenches  are  then 
caused  by  interference. 

We  have  in  the  past,  and  also  in  figure  3-1,  tried  to  follow  the  deepest  minimum 
and,  therefore,  found  some  strange  movements  of  the  dark  rings  and  sometimes 
extreme  difficulties  in  deciding  how  to  connect  the  minima.  The  suggested  inter¬ 
ference  of  the  two  systems  explains  our  difficulties.  It  shows  that  sometimes  the 
diffraction  effect  is  a  bit  larger  than  the  reflection  effect, and  the  overall  effect 
follows  the  solid  lines  in  figure  3-2;  at  other  times,  the  reverse  is  true.  In 
addition,  the  interference  can  also  lead  to  very  weak  minima  when  both  systems 
arrive  at  opposite  phase  but  equal  amplitude.  Furthermore,  the  meandering  of 
the  minima  in  the  backward  area  can  be  explained  as  the  interference  of  the  two 
systems  leading  to  a  deviation  from  the  straight  reflection  system.  If,  however, 
the  minima  follow  a  straight  line  with  no  meandering,  then  one  system  predom¬ 
inates  and  the  other  system  contributes  very  little. 

If  another  type  of  altitude  chart  were  used  in  which  the  angles  90*  to  180*  were 
distributed  according  to  another  formula  than  the  presently  used  one,  such  that 
8  =  180*  were  plotted  along  the  abscissa  forp  =  0,  then  the  diffraction  minima  in 
the  backward  area  would  follow  a  vertical  straight  line  similar  to  the  forward 
area.  Such  a  different  altitude  chart  might  be  useful  and  should  be  explored  to 
simplify  the  interpretation. 

3.  2  ALTITUDE  CHART  FOR  INTENSITY  FUNCTION  ^2 

The  pattern  on  this  altitude  chart  is  generally  similar  to  that  for  i^  although  less 
contour  lines  appear  (figure  3-4).  Hence,  most  of  the  statements  and  conclusions 
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made  for  apply.  The  description  given  for  n  =  1.  33  is  also  valid  and  need 
not  be  repeated.  Here  again,  figure  3 -5a  shows  the  minima  taken  from  figure 
3-4,  and  figure  3-5b  contains  the  two  intersecting  systems.  A  comparison  of 
figures  3-5a  and  3-5b  reveals  the  interaction  of  the  two  systems  as  discussed 
in  section  3.  2  for  i^. 

3.  3  POSITION  OF  DARK  RINGS 

Based  on  the  altitude  charts,  new  graphs  for  the  position  of  the  dark  rings  have 
been  constructed  and  are  shown  in  figures  3-6  and  3-7  up  to  a-  10,  These 
figures  supersede  the  curves  given  in  an  earlier  report  (Penndorf,  1961b)  and 
are  basically  similar  to  those  for  n  =  1.  33.  The  dashed  line  indicates  inflection 
points  only  and  no  real  minima.  In  addition,  the  two  systems  described  in  sec¬ 
tion  3.  1  are  added.  The  system  for  diffraction  minima  is  indicated  by  thin 
dashed  lines  going  from  the  lower  right-hand  side  to  the  upper  left-hand  side, 
and  the  system  for  reflection  minima  by  thin  dotted  lines  going  more  or  less 
perpendicular  to  the  first  system.  In  general,  the  actual  minima  follow  one  or 
the  other  system.  Meandering  indicates  that  both  systems  are  of  about  equal 
magnitude.  In  these  cases,  our  interpretation  of  the  continuous,  predominant 
minimum  may  not  always  be  correct,  and  a  more  detailed  investigation  based 
on  data  available  in  smaller  steps  than  A6  =  10°  may  change  our  interpretation, 
especially  around  9  =  90*.  Hence,  we  are  convinced  that  the  dividing  zone 
between  forward  and  backward  movement  is  not  of  great  importance  and  not 
fixed.  Thus,  it  is  believed  that  our  earlier  interpretation  that  the  dividing  zone 
is  situated  around  9  =  80°  to  90°  is  erroneous.  Data  analysis  for  other  refractive 
indices;  for  example,  n  =  1.  10,  1.  20,  1.  40,  and  1.  44,  support  our  present  view¬ 
point. 

Rowell  (private  communication)  has  completed  similar  studies  for  n  =  1.  486 
and  finds  a  much  more  complicated  pattern,  such  that  nearly  all  of  his  dark 
rings  move  to  larger  scattering  angles  with  increasing  size  parameter.  Hence, 
the  question  of  the  true  pattern  is  not  yet  solved.  Only  computations  in  smaller 
steps  of  9  can  prove  or  disprove  our  interpretation. 

3.  4  DIFFRACTION  MINIMA 

A  careful  investigation  of  the  diffraction  minima  can  also  lead  to  an  improvement 
of  the  scattering  diagrams  because  the  angular  position  and  the  depths  for  a 
particular  diffraction  minimum  changes  systematically  as  seen  on  the  altitude 
charts.  The  depth  of  the  valley  for  a  particular  a  value  is  defined  by 


namely,  the  ratio  of  at  the  lowest  point  of  the  valley  for  a  given  6  value  to  the 
value  i^  in  the  forward  direction  (6  =  0*).  The  parameter  D1  means  first  diffrac¬ 
tion  minimum.  An  analog  definition  applies  to  the  higher  orders  of  diffraction 
minima  (D2,  D3,  and  so  on)  as  well  as  for  i^* 

Data  of  such  i^^.  and  i-  have  been  computed  for  all  diffraction  minima  and 
plotted  in  a  coordinate’^system,  using  u  as  abscissa  and  log  i.  as  ordinate. 

For  the  first  diffracion  minimum,  data  up  to  a=  40  are  available;  for  the  rest, 
only  up  to  a=  10.  Since  all  diagrams  are  similar,  only  that  for  D1  is  shown  as 
an  example  in  figure  3-8  because  it  is  the  most  complete  one.  The  data  for 
os  10  to  o  =  19  are  less  certain  than  the  rest  because  ^8>  10*  in  this  size  range. 
Therefore,  they  are  left  out,  also  the  data  from  o=  2  to  cy=  5.  5  are  omitted 
because,  for  such  small  ^values,  the  pure  diffraction  pattern  is  not  yet  estab¬ 
lished. 

From  o=  5.  5  to  6.  5,  a  counterclockwise  motion  begins  followed  by  a  clockwise 
motion  to  o=  9.  5.  The  next  curve  begins  at  q=  19  in  a  counterclockwise  motion 
to  19.  6,  is  followed  by  a  clockwise  motion  to  (y=  23.  4,  and  so  forth.  The  minima 
for  u  >  3.  83  are  lower  than  those  for  u  <  3.  83.  The  lowest  values  of  i^^  (the 
trenches  on  the  altitude  chart)  occur  always  when  u  reaches  a  maximum  or 
minimum  for  the  counterclockwise  motion.  This  simple  system  of  clockwise 
and  counterclockwise  motion  has  been  proven  first  for  ns  1.  33  and  can  be  estab¬ 
lished  now  also  for  n  =  1.  5  and  to  higher  values  of  P.  It  has  been  suggested  in 
the  report  for  ns  1.  33  that  the  relation  given  in  section  3.  1  exists,  and  this  is 
now  proven.  For  all  n's,  the  minima  of  i^^  follow  the  relation 

Pt  =  2iT(y  ±  1/4) 

except  for  y  s  l,  and  some  slight  deviation  from  this  expression  may  be  real. 

For  the  other  diffraction  minima,  the  same  relation  should  hold,  but  more  data 
in  smaller  steps  of  9  are  needed  to  prove  it. 

3.  5  SCATTERING  DIAGRAMS 

To  visualize  the  numerical  results  and  to  show  the  correct  interpolation  of  the 
tabulated  data,  a  graphical  presentation  of  the  intensity  functions  i^,  i.i  and 
i-  +  12  has  been  prepared  in  the  form  of  scattering  diagrams  for  a=  0  to  10  steps 
of  A0f=  0,  5;  they  are  given  in  appendix  6.  1.  The  logarithmic  scale  chosen  was 
such  to  allow  reasonable  interpolation. 

The  method  of  construction  is  the  same  as  for  n  =  1.33  and,  therefore,  not 
repeated. 
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This  chart  represents  the  forward  area  for  9<  10  only.  The  empha¬ 
sis  is  placed  on  the  first  diffraction  minimum,  and  some  contours 
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The  position  of  the  dark  rings  is  indicated  by  solid  lines, 
resent  the  beginning  of  a  new  ring.  The  center  position  ( 
cated  by  full  circles. 

63-3082 
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Figure  3-8  RELATIVE  INTENSITY  FUNCTION  iij,  FOR  THE  FIRST  DIFFRACTION 
MINIMUM  AS  FUNCTION  OF  DIFFRACTION  PARAMETER  u  FOR  n  =  1.  5 

Minima  are  reached  when  u  approaches  extreme  values.  Numbers  on  curve  rep¬ 
resent  selected  o  values.  Data  are  based  on  altitude  charts  prepared  up  to  40. 
Heavy  dashed  line  represents  u  =  3.  83,  the  classical  value  for  the  first  diffrac¬ 
tion  minimum. 

63-3085 
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IV.  CONCLUSIONS 


Twenty  scattering  diagrams  of  the  intensity  functions  ii,  iz,  and  i^  +  iz  are 
given  in  the  form  of  an  atlas  for  n  =  1. 5  and  a  =  0. 5  (0.  5)  10.  They  are  based 
on  computations  in  steps  of  Aa  =  0.1  and  6  =  0  (1)  10  (10)  180*. 

Since  a  linear  interpolation  is  not  allowable  for  a  >  3,  interpolation  methods 
have  been  investigated.  The  altitude -chart  technique  allows  to  interpolate  the 
functions  graphically  for  any  set  of  a  -  8  values.  Furthermore,  the  position 
and  depth  of  the  minima  are  determined  correctly.  The  dark  rings  may  clearly 
be  divided  into  the  two  grovq>s:  those  due  to  diffracted  light,  and  those  due  to 
reflected  and  refracted  light.  For  the  first  group,  the  diffraction  parameter  u 
for  each  particular  dark  ring  meanders  around  the  classical  diffraction  param¬ 
eter  which  it  approaches  asymptotically  for  very  large  spheres.  For  the  second 
group,  the  diffraction  parameter  u  increases  proportional  top.  The  center  of 
the  deep  trenches  occur  for  the  first  group  in  steps  of  Apsa  V,  for  the  second 
group  in  steps  of  Ap  =  0.81.  The  two  systems  can  be  separated,  and  the  inter¬ 
action  of  the  two  systems  determines  the  direction  in  which  the  predominant 
minimum  moves.  This  interaction  also  explains  the  position  of  the  trenches 
and  the  areas  where  the  minima  suddenly  become  indistinguishable. 

The  scattering  diagrams  have  been  constructed  using  all  available  interpolation 
techniques  in  addition  to  the  computed  data  and  represent  the  most  reliable  data 
up  to  a  s  10.  They  allow  interpolation  for  any  set  of  a  -8  values. 
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INTENSITY  FUNCTION 


INTENSITY  FUNCTION 


6.  2  TABLE  OF  DARK-RING  PARAMETERS  FOR  INTENSITY  FUNCTIONS  ij 
AND  i2.  n  =  1.5  AND  a  =  0.5  (0.5)10 

Symbols  used  in  tables  6-1,  6-2,  and  6-3  are  as  follows; 


a 

P 

a 

u 

*lmin 

Mr 

D 

T 


Size  parameter 

Normalized  size  parameter 

Scattering  angle  in  degrees  (0*  =  forward) 

Diffraction  parameter 
Intensity  function  at  minimum 

Relative  intensity  function  for  particular  dark  ring 
Dark  ring  due  to  diffracted  light 

Dark  ring  due  to  transmitted  (reflected  and  refracted) light. 
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TABLE  6-1 


DATA  FOR  INTENSITY  FUNCTION  ij  AND  D1 


6 

u 

^Imin 

a=p 

8 

u 

llmin 

ilr{Dl} 

2.30 

3.6-3 

3.8-“* 

26.0 

7 

3.33 

1.5^ 

1.2-3 

wm 

2.50 

3.2"^ 

1.0-2 

26.5 

7 

3.16 

1.52 

1.1-3 

74 

2.88 

1.5:1 

1.9-^ 

27.0 

7 

3. 10 

6.o2 

4.3-3 

61 

3.05 

1.2-2 

27.5 

7 

3.18 

1.53 

1.1-2 

53 

3. 19 

3.5° 

1.3-2 

28.0 

7 

3.33 

4.o3 

2.1-2 

41 

2.89 

4.8° 

i.r2 

28.5 

7 

3.48 

3.o3 

1.4-2 

5.0 

35 

2.87 

I.5I 

2.5-2 

29.0 

8 

3.85 

1.53 

6.5-3 

5.5 

34 

3.07 

l.ll 

1.7-2 

29.5 

8 

4.05 

4.o2 

1.5-3 

6.0 

29 

2.91 

RIB 

2.5-2 

30.0 

8 

4.00 

1.52 

4.0-4 

6.5 

23 

2.54 

1.60 

2.3-2 

30.4 

7 

3.67 

4.o2 

1.4-3 

7.0 

19 

2.28 

IBl 

5.6-3 

31.0 

7 

3.68 

6.o2 

1.9-3 

7.5 

18 

2.32 

BBI 

5.7-^ 

31.4 

7 

3.62 

6.o2 

2.0-3 

8.0 

20 

2.73 

4.41 

5.4-2 

32.0 

6 

3.35 

3.52 

1.2-3 

8.5 

24 

3.45 

4.01 

3.1-2 

32.4 

6 

3.30 

3.o2 

1.0-3 

9.0 

28 

4.07 

7.8° 

3.7-^ 

33.0 

6 

3.30 

8.o2 

2.7-3 

9.5 

28 

4.60 

3.3O 

33.4 

6 

3.40 

2.5^ 

7.6-3 

mssM 

26 

4.37 

BOB 

9.5-* 

34.0 

6 

3.78 

5.0^ 

1.3-3 

34.4 

7 

3.92 

4.0^ 

1.0-2 

19.0 

11 

3.60 

3.52 

8.8-3 

35.0 

7 

4.02 

9.52 

2.0-3 

19.5 

10 

3.30 

4.0! 

1.0-3 

35.4 

7 

4.07 

6.02 

1.1-3 

20.0 

9 

3.00 

WESm 

2.4-3 

36.0 

7 

4.08 

5.0I 

8.7-4 

20.5 

9 

3.12 

2.52 

6.2-3 

36.4 

6 

4.05 

8.0I 

1.4-4 

21.0 

9 

3.21 

ElESa 

1.7-^ 

37.0 

6 

3.87 

1.52 

2.4-4 

21.5 

3.58 

2.6^ 

4.2-2 

37.4 

6 

3.73 

6.o2 

22.0 

3.94 

1.1^ 

1.5-2 

38.0 

5 

3.55 

9.o2 

1.4-3 

22.5 

4.18 

4.0^ 

4.4-3 

38.4 

5 

3.50 

5.o2 

8.2-4 

23.0 

4.23 

8.0I 

8.0-4 

39.0 

5 

3.53 

7.9^ 

1.3-3 

23.5 

3.92 

9.0I 

8.2-4 

39.4 

5 

3.57 

1.53 

2.5-3 

24.0 

3.96 

1.0^ 

7.7-4 

40.0 

5 

3.70 

5.o3 

7.5-3 

24.5 

9 

3.80 

1.62 

1.3-3 

25.0 

8 

3.65 

9.0I 

KKSB 

25.5 

8 

3.40 

4.0I 

3.1-3 
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TABLE  6-2 


DATA  FOR  INTENSITY  FUNCTION  ii.  DIFFRACTION 
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TABLE  6-3 


DATA  FOR  INTENSITY  FUNCTION  REFLECTION  SYSTEM 
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